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Spin Excitation in Nano-Graphite Ribbons with Zigzag Edges 
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Spin excitation in a nano-graphite ribbon with zigzag edges is investigated theoretically. Due 
to the strongly localized nature of the states near Fermi energy, the effective Hamiltonian for 
the low energy physics is given by Heisenberg Hamiltonian with the nearest neighbor exchange 
coupling. The action corresponding to the effective Hamiltonian is mapped to that of the 0(3) 
nonlinear sigma model. It is shown that the spin excitation has a gap when the number of the 
zigzag lines is even, whereas the excitation becomes gapless in case of the odd number of the 
zigzag lines. 
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Recently, the graphite-based one-dimensional materi- 
als with nano-meter sizes have been attracting much at- 
tention in both the fundamental sciences and the appli- 
cation sides. One of the most remarkable characteris- 
tics in these materials is that the electronic properties 
depend strongly on their geometrical structure. For ex- 
ample, the carbon nanotubes, which are made by rolling 
up a graphite sheet, can be either a metal or a semi- 
conductor depending on the chiral vector specifying the 
way of wrapping. 1-3 ) 

A nano-graphite ribbon (NGR) is the nano-meter size 
graphite fragment. It has been clarified that the edge 
regions play important roles for the electronic structure 
in this material. The calculation based on the tight- 
binding model revealed that a NGR with zigzag-shaped 
edges, which is called as a zigzag NGR in the following 
and shown in Fig.f, has strongly localized one-particle 
states around the edges. 4 ~ 6 ' The same conclusion is ob- 
tained also by the first-principle calculation. 7 ' Since the 
localized states appear near the Fermi energy, it is ex- 
pected that the mutual interaction between electrons af- 
fects the electronic states more strongly compared to the 
ordinary systems with extended states near the Fermi 
energy. Then it may be necessary for clarifying the elec- 
tronic properties to utilize the different approach from 
the method for the ordinary system. 

The effects of the mutual interaction on the zigzag 
NGR have been investigated by applying the mean-field 
theory to the Hubbard model with the on-site repulsive 
interaction and the hopping energy between the near- 
est neighbor atoms. 4,8 ' It has been shown that the in- 
finitesimal on-site repulsion induces spontaneous mag- 
netic moments; the conclusion is different from that of 
two-dimensional graphite sheet, where the finite value is 
needed for appearance of the magnetic order. In addi- 
tion, the large magnetic moments appear at the edges ( 
I A and NB sites in Fig.f ). Based on this spin struc- 
ture, the effective spin ladder model is proposed and the 
existence of the spin gap is claimed. 8 ' However, since it 



is considered that, in one-dimensional systems, the ef- 
fects of the mutual interaction are not sufficiently taken 
into account in the mean-field treatment, the conclusion 
obtained by it has to be reexamined. In the present 
work, we investigate the spin excitation of the zigzag 
NGR with N (number of the zigzag lines from one side to 
the other, see Fig.f) based on the renormalization group 
analysis and by mapping the 0(3) nonlinear sigma model 
(NLSM). Quite recently, Hikihara et al. investigated the 
ground state and the excitation around it of the zigzag 
NGR by both the bosonization theory and the density 
matrix renormalization group method. 9 ' The bosoniza- 
tion theory insists that the ground state is a spin-singlet 
Mott insulator with finite charge and spin gap. The nu- 
merical calculation supports the conclusion in case of 
N = 2. We compare our result with the above conclu- 
sions obtained by the previous studies. 




Fig. 1. The structure of the zigzag NGR with TV = odd. Here the 
close (open) circle shows the A (B) sublattice and the rectangle 
with the dashed-dotted line indicates the unit cell. 
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As a model of the interacting electrons in the zigzag 
NGR, we consider the Hubbard Hamiltonian used in refs. 
4, 8 and 9. At first we derive the effective Hamiltonian 
describing the low energy excitations. As was already 
discussed above, the states near Kp = n/a (a : lattice 
spacing) have the strongly localized nature around the 
edges, and then the energy dispersion near Kp are given 
as E(K) ~ ±t(Ka - ir) N as was shown in Fig.2. 6 ) When 
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Fig. 2. The band structure of the zigzag NGR with N = 2 (a) and 
N = 3 (b). Here t is the hopping energy between the nearest 
neighbor sites and a is the lattice spacing. The thin dashed 
curve express the Fermi energy and the asymptotic behaviors 
near Fermi energy, E(K)/t ~ ±(Ka — tt) n are written by the 
thick dotted curves. 



the effective low energy theory is constructed with tak- 
ing account of the states near Kp, the interaction terms 
including the four fermion operators have the scaling di- 
mension, N — l, at the tree level as was already discussed 
in ref.9. Therefore, the interaction is always relevant for 
N > 2. In this case, one electron is localized at the one 
site and the system behaves as Mott insulator. Thus, 
the charge degree of freedom is frozen, i.e., the charge 
excitation becomes gapped, and the low energy physics 
is determined by the spin degree of freedom. In such 
a case, we should choose an approach from the atomic 
limit. Then, the effective Hamiltonian is given by the 
Heisenberg model with the nearest neighbor exchange 
coupling J, which is written as follows, 



N L 



H = jJ2{ h odd(l) + h CVCD (l)} , 
i=i 

WO = E {SiA(i)s iB (i) + s iB (i - i)s iA (i) 



(i) 



i:odd 



+S iA (l)S i - 1B (!)}, (2) 

/Wn(0= Yl {SiB(l)S iA (l)+S iA (l-l)S iB (l) 



+S iA {l)Si- 1B {l)}, 



(3) 



where N L — L/a with L being the length of the NGR, 
S a/b{1) = S N+1A/B (l) = and S 2 iA/B {l) = S(S + 1) 
with S = 1/2. Though the magnitude of the spin, S, 
is a one-half, we do not restrict it in the following for 
convenience. 

The action A = Awz + Ao corresponding to the above 
Hamiltonian is expressed by the coherent state path in- 
tegral as, 10 ) 



N 



A W z = iSy^y^{w[n iA (l)] +w[n iB (l)}}, 
JS 2 



(4) 



i=i 
?2 rP 



An = 



I drJ2[Y{(n iA (l)+n iB (l)Y 

J ° l i:odd 



+ {n iB (l - 1) + n lA (l)) 2 + (n tA (l) + n^ 1B (l)) 2 } 
+ J2 {(n lB (l)+n lA (l)) 2 + {n lA (l-l)+n tB (l)f 



+ (n lA (l) + n l - 1B (l)f } 



(5) 



where (3=1/ (k B T) (T : temperature) and the constraint 
n ?A/s(0 = 1 i s imposed. The quantity, w[n iA / B {l)] is 
the solid angle which n iA / B (l) forms in the period < 

T<(3. 

In the classical limit, the Hamiltonian (l)-(3) has the 
ground state with Si A (l) = Se z and Si B (l) = —Se z 
where e z is the unit vector along the ^-direction. Here we 
discuss the fluctuation around the classical solution by 
expanding as Si A (l) = Se z +s iA (l) and 6^(7) = —Se z + 
Si B (l). The equations of motion of the fluctuation in the 
linearized approximation are given as follows, 



sf B Q) 



NB 



(0 



-US{2s+ A (l) + s+ B (l) + s+ B (l - 1)} , (6) 
-US{3st A (l) + st B (l) + sf_ 1B (l) 
+s+ B (l + (-lf)}, forz = 2~iV 
US{3 S + B (l) + S + A (l) + 4 +1A (l) 
+4 A ( l -(- 1 Y)}, forz = l~iV 
US {2s+ B (l) + s+ A (l) + 4a(1 (-1)")} , 

(9) 



(7) 
1, (8) 



where 
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iA/B 



(I) 



(I). 



By in- 



troducing the Fourier transformation, s^ A (l) = 



{Nl)- 1/2 T, 

{Nl)- 1/2 Y. 



j\kU+(-iy /A)a-wi\ „+ 
°i/ 

„i[fc(£-(-l)74)a-wt] + 



(fc) and s+ B {l) 



sJ B (k), the eqs.(6)-(9) 



are rewritten as, 



UJS 



l/i 



"4a 



(k) 
(k) 



2s+ 



(fc) + 2cos(fca/2)s+ B (fc), (10) 
3st A (k) + s+_ 1B (k) + 2cos(ka/2)s+ B (k), 

for i = 2 - N, (II) 



us+ B (k) = - {3 S + B (k) + , 



,+ 



u)s+ B (k) 



i __ lA (k) + 2cos(ka/2)s+ A (k)}, 

(12) 
(13) 



for i = 1 ~ N - 1, 
{2s+ B (k) + 2cos{ka/2)s+ A (k)} 



The above equations have the two Goldstone modes 
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whose frequencies are given by u> — uj/JS = 
±2|sin(fca/2)| - ±|fca| = uj±{k)/JS (see Fig.3). Since 
the eigenvectors for the modes are obtained up to the 
fc-linear as 

st A {k) = -{l± N ~l + l \ka\), (14) 

sf B (k) = l±^\ka\, (15) 
the fluctuation, s^ A / B (l), can be written as follows, 




Fig. 3. The dispersion relation of the spin wave of the N = 2 (a) 
and N = 3 (b) zigzag NGR. Here the only the positive values 
are shown. 
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\ka\)(3(k)}, 



\ka\) 
(16) 

(17) 



where a(k) = —\uj + (k)a(k) and $(k) = — \uj_(k)j3(k). 
When we introduce the spatially slowly varying functions 
made from the Goldstone modes, 



F(xi) 



1 \^ c ifc(i-l/4)q 



x {-(1 + \ka\/+)a(k) - (1 - \ka\/4)0(k)} > (18) 
G(xi) = ^ E c ifc( '- 1/4)a {-\ka\/2(a(k) - 0(h))} , 



(19) 



with xi = la, eqs.(16) and (17) are written as follows, for 



i =odd, 

st A {l)=F{x l ) + {N-i + l/2)G{x l ), (20) 
4(0 = ~F{xi + a/2) -(N-i- l/2)G(xi + a/2), 

(21) 

and for i =cvcn, 



Va(0 = F i x i + °/ 2 ) + ( N - 1 + 1 / 2 )G(xi + a/2), 



at B (l) = -F(x l )-(N-i-l/2)G{xi). 



(22) 
(23) 



Next, we try to map the action given by eqs.(4) and 
(5) to that of the 0(3) NLSM. Eqs. (20)-(23) suggest 
the following ansatz, for i=odd, 



n iA (l) = n(x l ) + (N-i + l/2)-l(x l ), 



(24) 



ni B (0 = -n(a;i+a/2)-(iV- 



and for i =evcn, 



l/2)-l(x l+ a/2), 
(25) 



n lA {l) = n( Xl + a/2) + (N-i + l/2)-I(x, + a/2), 



n iB (l) = -fl(xi) - {N 



i-l/2)-l( Xl ). 



(26) 
(27) 



The constraint In^/sl = 1 leads to |fi| = 1 and fll = 
to the first order of a/S. By substituting eqs.(24)-(27) 
into eq.(4), we obtain in the continuum limit, 



A 



wz 



WQ 



(28) 



where = 2irS for N —odd and 9 = for N =even. 
The quantity Q is the winding number of the Euclidean 
space configuration {fl(x,r)} and written as, 

Q = -!- / drdx(d T n x d x £l)n. (29) 

47T J 

On the other hand, the quantity Aq is expressed as 

A Q = J drdx {l(d x nf + |/ 2 }. (30) 

By integrating the field I in eqs. (28) and (30), we obtain 
the following effective action of the NLSM as 

1 



2ff 



d'x{d^y 



(31) 



where X\ = x, x 2 = vt with v = JSa and g = 
2 / (NS) . The above effective action is the same as that 
for spin ladder systems 11-13 ) and the excitation is gap- 
less (gapped) for N — odd (even) in case of S = 1/2. 
The conclusion is different from the result based on the 
mean-field theory and that by the bosonization theory 
where there is the spin gap irrespectively of the number, 
N. The discrepancy is considered to be due to the fol- 
lowing facts. The effects of the interaction are difficult 
to be taken into account sufficiently in the mean-field 
theory starting from the band picture and the consistent 
mapping from the Fermionic theory to the bosonic one is 
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not done in the present Case without the linear dispersion Brezin and J. Zinn-Justin (Elsevier Sciences Publishers B.V., 

across the Fermi energy. In ref. 9, the spin excitation 1989) 

has been calculated for N — 3 numerically. It seems to 

us that the result indicates that the excitation for N = 3 

is gapless, which is identical with our conclusion, though 

the authors do not insist it because of lack of the enough 

system size for extrapolation. 

We note that the velocity of the excitation v is iden- 
tical with that obtained by the spin wave analysis. In 
addition, in case of N = 1, the present results are iden- 
tical with that with one-dimensional spin chains 14 ) with 
a — > 2a because the lattice spacing becomes a half. 

The coupling constant, g, for TV = 2 obtained here is 
smaller than that for the two-leg spin ladders with the 
isotropic coupling. 12 ' 13 ) Then the spin gap of the present 
system is expected to be smaller than that of the spin 
ladders. This fact is considered to be due to the lack of 
the exchange coupling between eg., 1A and 2B. In de- 
riving the effective action, eq.(31), the modes with high 
frequencies have been neglected. We have derived the ef- 
fective action with taking account of both the Goldstonc 
modes and the high frequency modes in case of N = 2. 
The conclusion is the same as eq.(31) with N = 2. Thus, 
it is expected that the high frequency modes do not play 
an important role for the analysis presented here. 

In summary, we investigated the spin excitation of the 
zigzag NGR with the width N. We concluded that the 
spin excitation is gapless for N = odd, whereas the ex- 
citation has a gap in case of N = even. 
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